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Abstract 

It is derived that the two-dimensional Ising model with alternat- 
ing/random interactions and with periodic/free boundary conditions is 
equivalent to the ground state of the one-dimensional alternating/random 
XY model with the corresponding periodic/free boundary conditions. 
This provides an exact equivalence between a random rectangular Ising 
model, in which the Griffiths-McCoy phase appears, and a random XY 
chain. 
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Let us consider the two-dimensional rectangular Ising model 
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where {crfj}(k 



x, y, z) are the Pauli operators satisfying <J^ +M j = afj 
of- . The partition function of the model ([!]) is obtained 



and °ij+N 

from the maximum eigenvalue of the transfer matrix 
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V 2 = exp[^X>?aJ +1 ], K) = j]/k B T (l = v,h). 



Here {cr^} are again the Pauli operators satisfying cr^ +M = <jj, and 
KV* is the dual interaction of K] defined by tanhiCJ* = exp(-2i^J). 

Suzuki derivedpQ [2] that the transfer matrix with the uniform 
interactions KJ = K v and A'j 1 = K h commute with the Hamiltonian 
of the one-dimensional quantum XY model 
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with the uniform interactions — J k (k — x, y) and the uniform 
external fields Zfj = if. He showed that two operators V and ffxY 
commute when the coupling parameters satisfy the condition that 

J y /J x =t&nh 2 K m , (4) 
HH/J X = 2tanhfT ; *coth2fi:' i . (5) 

Hence V and ifxY can be diagonalized with a common basis set of 
eigenvectors. In particular, the eigenstate for the maximum eigenvalue 
of V coincides with the ground state of ifxY , and hence the thermody- 
namic properties, such as correlation functions and critical singularities 
of the two-dimensional Ising model are related to those obtained in the 
ground state of the one-dimensional quantum XY model. 

In this letter, it is derived that this equivalence can be extended to 
the models with general coupling parameters including alternating and 
random interactions and external fields with periodic or free boundary 
conditions. The condition 



[V,if XY ] =0, 



(6) 
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is equivalent to 

V 2 V 1 1/2 H XY V 1 - 1/2 V 2 - 1 = V- 1I2 H^V1 12 . (7) 
One obtains 



N 
3=1 



+ JJ(- Cj ofoJ +1 + C+oJoJ +1 - Sj iofo? +1 - S+iof ojf +1 )] 

/^X, (8) 
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where 
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The coefficient of the term <Jj<Jj +1 vanishes if we assume 

<<', = ^4- ( 9 ) 

Then it is straightforward to derive that ([7]) is satisfied if we assume 



fj,Hj 
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where = JJ — sjjj, and 
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-[cosh(2i<'j l _ 1 + 2^) ± cosh(2^ l _ 1 - 21$)] 



_ ± = l smhQK^ + 2itf ) sinhggLi - 2K>J) 

and the symbols [H — ]j and [ — h]j denote 

[±,T] i = 2^_ 1 s J ± +2^sJ. 

From the equation (J9j> , JJ and Jj can be written as 

Jf = (cosh K]*) (cosh ^i)^-, 

JJ = (sinhi^Xsmhi^r,, (12) 
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Therefore, V and Hxy commute if (|T0[) and (1121) arc satisfied. 

There are solutions for the equations © and (|10p in the following 

cases: 

1-1) Uniform interactions: The interactions of the Ising model are 
uniform, ifj = if v and ifj 1 = K h . In this case, cf — , sf = s, 
2^ = c ± , and = s ± , where 

c ± = - [cosh 2if u * ± 1] , s = ~ sinh 2iT'* , 



i[cosh4if"±l], ^ = l[ Sinh4 f ±1]. 
2 L J 2 L 4if h 1 



Assuming the uniform interactions Jj = J k (k = x, y) and the uniform 
external fields Hj = H, the equations in (fT2"| result in and four 
equations in (ITU)) reduce to one identical condition 

which results in ([5]). 

1-2) Alternating-sign ifj 1 : Let us consider the case in which the 
vertical interactions of the Ising model are uniform, Kj — if 11 , and 
the horizontal interactions alternate in sign as ifj 1 = (— l) J if™. In 
this case, = c ± , = s, = icF* 1 , and = ±s ± . Assuming 



JJ = (-iyj h (k = x,y), i.e. tj = (-1) j t, we again obtain (fT3|) . 

1-3) Alternating if J: Let us consider the case in which the vertical 
interactions of the Ising model alternate as ifj = K\ for j = odd and 
KJ = if B for j = even, and the horizontal interactions are uniform, 
ifj = K h . In this case cf = c^ B , sf = s^ B for j = odd and = sJ B 
for j = even, cf — , and = s 1 * 1 , where c^ B and s^b are obtained 
from (fTT|) substituting ifj = if£ and ifj+i = K#. Assuming Tj = r, 
the equations in (|10p reduce to one condition that results in alternating 
external fields 

1 sinh 4/^ 

nHj = 2 (- sinh2ifr) r f . 

' J V 2 jy cosh4if h -l 

In this case, if, are alternating for j being even and odd, but JJ and 
JJ are uniform, i.e. JJ — J x = (cosh i^*) (cosh if g*)r and Jj = J* = 
(sinh ifX*) (sinh ifg*) r. We assume that if^ and ifg have the same 
sign so that the couplings of the XY model remain real (note that the 
dual interaction of -K v is (-K v )* = K v * ± ni/2). 

Two cases 1-2) and 1-3) are not exclusive. The Ising model with 
alternating-sign if j 1 and alternating if J has its equivalent XY chain. 

1-4) Random-sign ifj 1 : Let us consider the case where the horizon- 
tal interactions of the Ising model have the same absolute value but 
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random signs, K^ = £jK h , where ej = ±1. In this case we obtain 
= c*, = s, and also obtain that — and = for 
Cj-i = ej, and = ±c ± and = ±s ± for ej_i = — e,-. Assuming 
the sign of the interactions as Jj — ejj k (k = x,y), i.e. tj = CjT, we 



again obtain the solution ([13)) 

1-5) Arbitrary KJ: Let us consider the case where the vertical 
interactions of the Ising model {KJ} have arbitrary strength with the 

same sign. A solution exists provided that JJ_\ = J^ for each j. In 
this case, Tj — r in (TT2"j) and the external fields are obtained as 

sinh 4ft'' 1 i 

liHi = 2s,r , ' -, = -sinh2!n*. (14) 

p J J cosh^R^-l' J 2 J v ; 

The cases 1-2) and 1-3) are special cases of 1-4) and 1-5), respec- 
tively. Again two conditions 1-4) and 1-5) are not exclusive. Assuming 
1-4), Kj- — €jK h and Jj = tjJ k (k = x,y), a solution for 1-5) exists 

when JjL\ — Jt ■ The external fields are obtained from flT4"|) . 

For all of these cases, from 1-1) to 1-5), a solution exists 2-1) with 
the periodic boundary condition, and 2-2) with the free boundary con- 
dition. Assuming that K% = 0, we obtain = cos\v2K h , = 0, 
sjf = sinh 2K h /2K h and = 0. Assuming a finite XY chain tn = 0, 
and hence = 0, the condition ((9|) for j — N is satisfied with = 0, 
and (fit)]) for j = N have a solution 

sinh 2K h 

» H " = J N-X cosh2Kh _ r 

Some of these systems l-l)-l-5), with 2-1) or 2-2) have been intro- 
duced in connection with experiments and theoretical interests. 

Fisher and Ferdinand calculated]!] the free energies of two-dimensional 
Ising models on various lattices with various boundary conditions. 
Their motivation was to explain the shifts and rounding of the spe- 
cific heat maximum observed in experiments, in terms of the effects 
of microcrystalline structure. They introduced a grain boundary in 
which a vertical ladder of j horizontal bonds has modified interactions 
£J. Their case with £ = and j — > oo corresponds to our case 1-1) 
with 2-2). 

The two-dimensional Ising model with the periodic boundary con- 
dition in the horizontal direction and the free boundary conditions in 
the vertical direction was considered in a study by Schultz et al.[3], 
in which they diagonalized the model as a many-fermion system and 
obtained the free energy. Abraham solved[5] the eigenvalue problem 
for the transfer matrix of the rectangular Ising model with the periodic 
boundary condition in the vertical direction and the free boundary con- 
ditions in the horizontal direction, and investigated the detailed band 
structure. This model corresponds to the case 1-1) with 2-2). 
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The Ising model on a cylinder with the boundary condition that the 
spins on one edge are fixed to be up and the spins on the other edge are 
fixed to be down[B], and that with another boundary condition that the 
spins on the two edges are fixed to be up [TJ, were considered in order 
to shed light on the surface tension problem between two oppositely 
magnetized phases. 

McCoy and Wu[8] considered the rectangular Ising model with uni- 
form horizontal interactions, and vertical interactions that vary ran- 
domly from row to row, with the periodic boundary condition in the 
horizontal direction, and with the free boundary conditions in the ver- 
tical direction. They attempted to explain experiments that showed 
that the specific heat is a smooth function at T c . They suspected that 
this effect was due to the presence of random impurities. This model 
coincides with our case 1-4) with 2-2) if the distribution of the ran- 
dom interactions is assumed to be random ± J, though they assumed 
a narrow power law distribution function in their calculations. Their 
distribution function can be introduced in our case 1-5) with 2-2), and 
it is not unreasonable to expect similar behavior in our system. They 
obtained the result that the specific heat is not divergent and infinitely 
differentiable at and near T c though it is not analytic. It was also de- 
rived in another paper [5] that the usual parametrization in terms of 
critical exponents does not describe the critical behavior of this model. 

Shanker and Murthy[10][Tl] considered a model in which the verti- 
cal interactions are fixed and ferromagnetic, and the horizontal interac- 
tions are equal in each row and vary randomly in sign and in magnitude. 
Their model generally includes frustration. They assumed the periodic 
boundary conditions in two directions. This model corresponds to our 
case 1-5) with 2-1), provided that the horizontal interactions are ran- 
dom but have the same sign. They mapped the problem to a collection 
of one-dimensional random field problems, and identified three phase 
transitions. 

The commutation relation ^ certifies that V and Hxy can be 
diagonalized simultaneously. Moreover, it can be derived that the 
eigenstate for the maximum eigenvalue of V coincides with the ground 
state of Hxy, and hence the thermodynamic properties of the two- 
dimensional Ising models and the ground state properties of the one- 
dimensional XY models are related to each other. 

In order to show this coincidence, let us consider the general proper- 
ties of V and Hxy- The matrix Ylj=i ^j* a j ^ s diagonal, and hence all 
the elements of exp[^^ =1 ifJVJ] are non-negative. When Kj > 0, the 
elements of J2jLi Kj'&j&j+i are non-negative and hence the elements 

of exp[J2j =1 Kj<j'j<jj +1 ] are also non-negative. Therefore, in this case, 
all the matrix elements of V are non-negative. Because the Ising inter- 
actions are written as crjcrj +1 = sfst +1 +sjsj +1 +stsj +1 +sjs'j +1 , the 
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matrix V becomes block-diagonal and contains two irreducible block 
elements, one is which operates on the bases states with even 
number of up spins and the other is which operates on the bases 
states with odd number of up spins: 



V 



o 



o 

02) 



The block element is irreducible because (yW) m is irreducible for 
sufficiently large m. (This comes from the fact that (yW) m cannot be 
irreducible if is reducible.) Similarly is also irreducible. Let 



$j and be the eigenstates of V for the maximum eigenvalue of 
and V^ 2 \ respectively. These are expressed as 



,(2) 









Pi 



where (jyp and </>| 1 ' ) are vectors with 2 N /2 elements, and is the zero 
vector with 2 N jl elements. From the Perron-Frobenius theorem, the 



1 i ("2") 

elements of 4^' and 4>\ are non- negative. Because the two block 
elements and V ^ are irreducible, and (j)^ are non-degenerate 
eigenstates of and V^ 2 \ respectively, and the elements of cf>f' and 
(f>i are all positive. 

Next let us consider the Hamiltonian Hxy with Jj > (k = x, y). 
When we consider a matrix -ffxY+c/, where / is the unit matrix and c 
is a positive and sufficiently large constant. The eigenstates of — Hxy + 
cl are also the eigenstates of Hxy- In particular, the ground state of 
the Hamiltonian Hxy is the eigenstate for the maximum eigenvalue of 
—Hxy + cl. Because the XY interactions are written as 



JX X X 

J J 3 a 3 + l 



jV a V a V - 



the Hamiltonian Hxy is block-diagonalized into two blocks with even 
and odd number of up spins. From (|12l) . the interaction constants 
satisfy JJ > Jj, and hence all the matrix elements of —Hxy + cl arc 

non-negative and the two block elements are irreducible. Let $xy ano - 

(2) 

$xy b e t ne eigenstate of Hxy for the lowest energy eigenvalue of even 
and odd block element, respectively. These are expressed as 



(i) 

XY 



'XY 





(2) 
XY 





6 (2) 
9 XY 



It follows that the elements of 



^XY 



and 



4 (2) 

^XY 



are all positive. 
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The matrices V and Hxy are both Hermitian, and hence their 
common eigenstates form (or can be arranged to form) an orthogonal 
basis set. Therefore, any two eigenstates of V and Hxy are orthogonal 
or identical, or belong to a degenerate eigenspace. Because all the 
coefficients of <^>j and fixy are all positive, $j and ^xy carm °t be 
orthogonal. Considering the fact that (f>^p and 4>^ Y are non-degenerate 
eigenstates of the corresponding block elements, respectively, it follows 
that $j = except overall constant. Similarly, we obtain •^j 2 '' = 
<i> X Y except overall constant. 

In the cases with ifj 1 < and Jj < (fe = x, y) for some j, one can 
introduce a rotation of tt around the z axis in the spin space, and the 
problem is mapped to one with > and jj° > (k = x, y), which 
yields the same result. 

The remaining problem is to identify the eigenstate for the maxi- 
mum eigenvalue of Hi and the ground state of Hxy- For the periodic 
case 1-1) with 2-1), it has been shown in [5] through a direct diago- 
nalization by fermion operators that the eigenstate for the maximum 
eigenvalue of V is $j 1 ' > , and that the ground state of Hxy is 3?xy- 

In the case of 1-1) with 2-2), the XY model has free ends. The 
Hamiltonian Hxy is separated into two sectors, one with even number 
of fermions ( i.e. the block-element with even number of up spins), 
and the other with odd number of fermions ( i.e. the block-element 
with odd number of up spins). It is most clearly seen in the calculation 
by Niemeijer[12] that the energy eigenvalues and the eigenvectors are 
obtained from an eigenvalue problem for a symmetric matrix. Hence 
the XY model with free boundary conditions can be solved, and the 
ground state is unique and found in the even sector. 

Schultz et al. [3] diagonalized the two-dimensional Ising model with 
free boundary conditions in one direction. The maximum eigenvalue 
of the transfer matrix was found in the sector with even number of 
fermions ( i.e. in the block-element with even number of up spins), 
and was non-degenerate above T c . Below T c , the maximum eigenvalue 
of the odd sector becomes degenerate with that of the even sector in 
the thermodynamic limit. 

Consequently, we again find, for the case 1-1) with 2-2), that the 

eigenstate for the maximum eigenvalue of V is ^j 1 '', and the ground 

state of -ffxY is $xy' where = ®xy except overall constant. 

The cases with Kj < and J* < (k — x, y) for some j can be 
mapped to the cases with K 1 - > and Jj > 0(k = x, y) by a rotation of 
7r around the z axis in the spin space. Because the energy eigenvalues 
are invariant by the rotation, the results for 1-1) with 2-1) or 2-2) are 
still valid in the cases 1-2) and 1-4) with 2-1) or 2-2). 

Finally, let us consider the cases 1-3) and 1-5) with 2-1) or 2-2). 
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The system with K j 1 > can be mapped to the system with < 0, 
and the rectangular Ising model with KJ < and K^ < for all j can 
be mapped to that with KJ > and K!f > 0. It is easy to check that 
the corresponding XY models are also mapped each other by a gauge 
transformation. Therefore, we can assume that KJ > and Kj > 
for all j. Then the elements of V\ and V2 are all non-negative because 
KJ* are real and Kj 1 are positive. Without loss of generality, one can 
assume that Tj > in (TT2]) and then the interactions JJ and Jj satisfy 
Jj — > and Jj + JV > 0. All the non-zero elements of V and 
—Hxy + cl remain non-zero as long as KJ and Kj remain real, finite 
and non-zero. 

From the Perron-Frobenius theorem, the maximum eigenvalue of 
each block element is non-degenerate. When we assume that the eigen- 
values are continuous in terms of the coupling constants, this leads to 
the result that level crossings do not occur in each block element as 
long as all the non-zero matrix elements remain non-zero. Therefore 

0j = is the eigenstate for the maximum eigenvalue of the even 

(2) (2) 

block element, and (j>\ — (fr^y- is that of the odd block element, for 
real, finite and non-zero KJ and Kj. The state $^ = &xy ^ s ^ ne 
eigenstate for the maximum eigenvalue of V as well as the ground 
state of -Hxy, at least in a neighborhood of the uniform point 1-1). 

The correlation functions of the two-dimensional Ising model and 
those in the ground state ^>x Y of the one-dimensional XY model are 
related to each other. Let f(a) be a product of the Pauli operators aj. 
When one take the limit M — » 00 in the equation 

Tr/(q)(yiy 2 ) M _ Trf(a)Vi /2 V M V- 1/2 
Tr(l/il/ 2 ) M ~ TrV M 

one obtains the relation of the expectation values as 

</M>I = <V 1/2 / W 1/2 >XY, 

where (/(<x)}i is the expectation value of a spin product on the same 
row of the Ising model. In particular, the two-spin correlation functions 
satisfy the relation 

<a*af fc )! = cosh jq cosh K v k {a?a%) XY - sinh^J sinh!^>J^>xY, 

where we used the fact that all the elements of ^xy are rea ^ an d hence 
(ofoDxY = for j + k. 

As shown in |2] , the critical point T c of the two-dimensional Ising 
model corresponds to the critical field H c of the one-dimensional XY 
model, and two quivalent models show the same critical singularities. 
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However, especially in the cases with random coupling parame- 
ters, complex critical behavior can occur depending on the distribution 
function for {KJ}. It is known that the Griffiths-McCoy phase, which 
was originally found in two-dimensional random Ising models. [5] 1 13] 
appears in the one-dimensional random transverse Ising model. [T4"] 
and one-dimensional XY model with random interactions and ran- 
dom fields [T5]. As mentioned above, it is also found[10][TI] that the 
Griffiths-McCoy phase appears in the two-dimensional Ising model 
with random horizontal interactions, which corresponds to our case 
1-5) with 1-1). The relations derived in this letter shows the exact 
equivalence of certain kind of two-dimensional random Ising models 
and the one-dimensional random XY models. 
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